The role of surface plasmons in the decay of image-potential states 

on silver surfaces 
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The combined effect of single-particle and collective surface excitations in the decay of image- 
potential states on Ag surfaces is investigated, and the origin of the long-standing discrepancy 
between experimental measurements and previous theoretical predictions for the lifetime of these 
states is elucidated. Although surface-plasmon excitation had been expected to reduce the image- 
state lifetime, we demonstrate that the subtle combination of the spatial variation of s-d polarization 
in Ag and the characteristic non-locality of many-electron interactions near the surface yields sur- 
prisingly long image-state lifetimes, in agreement with experiment. 
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The understanding of charge locahzation and electron 
dynamics at solid surfaces is of crucial importance in a 
large variety of physical and chemical phenomena, such 
as energy transfer, electronically induced adsorbate reac- 
tions, catalytic processes, epitaxy, and the development 
of new materials [|l[ ||, |^, ^. Accordingly, over the last 
two decades there has been a large effort to investigate 
the characteristic energies and linewidths of short-lived 
electronic excitations at metal surfaces ^, H, 0, ^ • 

A special class of electronic excitations at metal sur- 
faces are the so-called image-potential states, which are 
localized on the vacuum side of the crystal but can move 
almost freely in a plane parallel to the sur face |, 0. 
These surface states occur when an electron is trapped 
between its own image potential on the vacuum side and 
a projected band gap of available bulk electronic states 
on the crystal side of the surface. Because of their well- 
defined physical properties, image states are useful to in- 
vestigate the coupling of excited surface electronic states 
with the substrate in the time domain [|ll|, |l^ . 

Recently, femtosecond ultraviolet-pump and infrared- 
probe techniques were combined with two-photon pho- 
toemission (2PPE) to provide accurate measurements of 
the relaxation dynamics of image-potential states on the 
(100) surfaces of Cu and Ag @ |l|]; the lifetime of the 
n = 1 image state on these surfaces was found to be 40 ±6 
and 55 ± 5 fs, respectively, corresponding to linewidths 
of approximately 16.5 and 12 meV. Recent theoretical 
work based on many-body calculations of the electron 
self-energy, which revealed the main factors that deter- 
mine the coupling of image states to the substrate, pre- 
dicted image-state lifetimes that are in the case of Cu 
0, |l6|, and other metal surfaces (l|, |l|, in good 
agreement with experiment but failed to explain the mea- 
sured relaxation of image-state electrons in Ag [EOl . 



In this Letter, we investigate the decay of the Ag(lOO) 
71 = 1 image state by treating Auger-like single-particle 
and collective surface excitations on the same footing. A 
consistent treatment of these decay channels is particu- 
larly important because of the near-degeneracy of the im- 
age state energy with the Ag(lOO) surface plasmon. The 
surprising and novel result of our work is that, although 
the imaginary part of the image-state self-energy is en- 
hanccnd due to the presence of the plasmon decay chan- 
nel, interferences resulting from the non-locality of the 
self-energy lead to a smaller overall image-state broaden- 
ing, in agreement with experiment. 

While measured image-state lifetimes are found to be 
systematically longer in Ag than in Cu surfaces Q, 
previous theoretical estimates seemed to indicate that 
image-state lifetimes should be shorter in the case of 
Ag. Although image states on Ag surfaces have ener- 
gies that are slightly closer to the Fermi level than in 
the case of Cu, both a smaller density of sp valence elec- 
trons and a smaller polarizability of d bound electrons in 
Ag are expected to yield weaker screening and therefore 
shorter image-state lifetimes in this material, as occurs 
in the case of bulk excitations |21, 22 . Furthermore, 



image-state lifetimes in Ag were expected to be further 
reduced, due to the opening of the surface-plasmon exci- 
tation channel [p3|. As we argue below, this anticipated 
trend is reversed once the crucial nonlocality of the self- 
energy in the surface region is taken into account. 

A first-principles description of electron dynamics in 
the bulk of noble metals shows that deviations from elec- 
tron dynamics in a gas of sp electrons mainly originate 
in the participation of d electrons in the screening of 
electron-electron interactions [^, hence, we account 
for the presence of occupied d bands in Ag via a po- 
larizable medium giving rise to additional screening, as 
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suggested in Ref. |g6[ , and compute the hfetime from the 
quasiparticle self-energy. 

The inelastic lifetime of image-state electrons at a solid 
surface is known to be determined by electron-electron 
many-body interactions. On the energy-shell (i.e., ne- 
glecting the quasiparticle energy renormahzation) and 
assuming translational invariance in the plane of the 
surface, which is taken to be normal to the z axis, 
the decay rate (or lifetime broadening) of the excited 
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tained by solving the one-dimensional Schrodinger equa- 
tion modeled in Ref. [^. To derive v' we assume that 
the sp valence electrons are embedded in a polarizable 
background characterized by a z-dependent local dielec- 
tric function ed{z,uj) = ed{uj) for z < Zd and ed{z,uj) = 1 
for z > Zd, and find p3| 
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state 4>{z) e*''" ''!! with parallel momentum k|| and energy 
E = £-|-fcjj/(2m) (unless stated otherwise, we use atomic 

units, i.e., e^ = Ti = me = 1) is given by p] 



q\\ <^diz',i^) 
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= -2 dz dz' (f>*{z)lraJ:i, ,{z,z')(l){z'), (1) via e(w) = e,p(a;) -h ed(w) - 1 (see, e.g., Ref. Q) 
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where 'E,]s^^^ ,;{z,z') is the two-dimensional Fourier trans- 
form of the quasiparticle self-energy and m represents 
the effective mass accounting for the potential variation 
in the plane parallel to the surface. 

In the GW approximation one considers only the 
first-order term in a series expansion of the self-energy E 
in terms of the screened Coulomb interaction W. If one 
further replaces the Green function G by that of nonin- 
teracting electrons one finds 



where ad — [^di^^) ~ ll/I^dl^) + !]• <^di^) is related to 
the long-wavelength limit of the bulk dielectric function 

where 



Imi;k||,e(^:, z') 
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where (/)/(2:) e*''ii''"ii represent the available single-particle 
states with energy Ef — Sf + (k|| — q||)^/(2m/) {Ep < 
Ef < E), Lu — E — Ef is the energy transfer, and Ep 
is the Fermi energy. Short-range exchange-correlation ef- 
fects, which are absent in the GW self-energy of Eq. (||), 
can be included in the framework of the so-called GWT 
approximation [^. GWT calculations for the n = 1 
image-state lifetime on the (100) and (111) surfaces of 
Cu were reported in Ref. [p^ , showing that they are very 
close to GW calculations as long as the screened interac- 
tion W is treated within the random-phase approxima- 
tion (RPA). 

Here we report GW-KPA calculations of the electron 
self-energy and the lifetime broadening by combining 
the one-dimensional model potential used in previous 
work |l^, |l^, which accurately accounts for the dy- 
namics of sp valence electrons p9|] , with the model re- 
ported in Ref. p3t in which the occupied d bands are 
accounted for by the presence of a polarizable medium 
which extends up to a certain distance from the surface. 
Within this model, the screened interaction is of the form 
W = v' I (1— x°w')i where v' is the d-screened Coulomb in- 
teraction and ^ is the non-local dynamical susceptibility 
of noninteracting sp electrons. We emphasize that this 
approach provides a coherent treatment of single-particle 
and collective surface excitations. 

The single-particle wave functions entering Eqs. (|l|) 
and (^ , and those entering the expression for x° , are ob- 



esp(a;) corresponds to the long- wavelength limit of the 
bulk dielectric function of sp valence electrons. The co- 
ordinate Zd represents the position of the boundary up 
to which the polarizable medium is assumed to extend. 

For comparison, we note that in a spatially uniform 
system the polarizable medium extends to infinity (z^ — ^ 
oo) and the d-screened Coulomb interaction is simply 
v' = V /ed{^)i V being the bare Coulomb interaction. The 
exact high-density limit of Eq. (|l|) is then found to be, for 
electron energies near the Fermi level [E — Ep « Ep), 
= T~p/y/€d{u) 0), as suggested by Quinn [|l|, 
where tqp represents the formula first derived by Quinn 
and Ferrel for a free-electron gas |^ . 

The n = I image state on Ag(lOO) is located close to 
the center of the projected band gap, it has a binding 
energy of 0.53 eV, and e — Ep = 3.9 cV. Its probability- 
density has a maximum at 3.8 A outside the crystal edge 
{z = 0), which we choose to be located half a lattice spac- 
ing beyond the last atomic layer, and the penetration into 
the bulk crystal (z < 0) is found to be 5%, as in the case 
of Cu(lOO). The effective mass of the n = 1 image state 
on Ag(lOO) is known to be very close to the free-electron 
mass (m ^1), and for all single-particle energies enter- 
ing x° and Im[— E] we choose the effective mass m/ to 
increase from our computed value of 0.45 at the bottom 
of the gap to the free-electron mass at the bottom of 
the valence band. The local dielectric function ed{i^) is 
taken from bulk optical data The boundary of the 
polarizable medium is taken to be at = — 1-5 ao (ao is 
the Bohr radius) inside the crystal, which was previously 
found to best reproduce the anomalous dispersion of Ag 
surface plasmons [ p6| . 

The lifetime broadening of the n — 1 image state at 
the r point on Ag(lOO) originates in processes in which 
the image-state electron with energy E = Ep + 3.9 eV de- 
cays into an empty state with energy Ef above the Fermi 
level {Ep < Ef < E). These processes can be realized 
by transfering energy and momentum to an excitation of 
the medium, thereby creating either a single-particle or 
a collective surface excitation with energy oj = E — Ef 
(0 < w < 3.9 eV). Without d-electron screening the 
surface-plasmon energy is too large {tOsp = 6.5 eV) and 
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FIG. 1: Imaginary part of the n = 1 image-state self-energy 
Im[— Ek|| ,2(2, 2')]i versus z, in the vicinity of the (100) surface 
of Ag . The value of z (indicated by an open circle) is fixed at 
-10 (top panel), (middle panel), and 10 ao (bottom panel). 
Notice the different scales. Solid and dashed lines represent 
the results obtained with zd = — 1.5ao and with zd —00, 
respectively. Full circles represent the atomic positions. The 
geometrical electronic edge (z = 0) is chosen to be located half 
an interlayer spacing beyond the last atomic layer. Parallel 
momentum and energy of the quasiparticle are taken to be 
k|| = and e = _Bf + 3.9 eV. 

only single-particle decay takes place; the presence of d 
electrons reduces Wgp to ~ 3.7 eV so that relaxation via 
surface plasmons becomes feasible. 

Fig. shows the imaginary part of the self-energy 
Im[— I]k||,£(2;, 2;')] foi" l*|t = and e — Ep = 3.9 eV, corre- 
sponding to the n — 1 image state at the T point of the 
Ag(lOO) surface, as a function of z and for fixed values 
of z'. In the top panel z' is taken to be well inside the 
metal (z' — — lOoo), showing that Im[— E] has a maxi- 
mum at z = z', as in the case of a homogeneous electron 
gas. When z' is fixed at the crystal edge {z' = 0), as 
shown in the middle panel of Fig. Im[— S] is still found 
to be maximum at z = z', but the magnitude of this 
maximum now being enhanced because of the reduced 
electronic screening near the surface. The bottom panel 
of Fig. [1] corresponds to z' outside the metal (z' = 10 ao). 
In this case, the main peak of Im[— E] lags behind and 
remains localized at z rather than z = z', showing a 
highly non-local behaviour of the imaginary part of the 
self-energy in the presence of a metal surface. The largest 
values of Im[— E] are obtained for z, z' « 0, because of 
the high probability of single-particle and collective elec- 
tronic excitations in the surface region. 

The imaginary part of the image-state self-energy in 
the absence of d electrons is also represented in Fig. |lj as 
obtained with Zd — > —00 [or, equivalently, ed{uj) = 1]. 
Although d electrons giving rise to additional screen- 
ing should reduce both W and Im[— E], this reduction 
is outweighted by the opening of the surface-plasmon ex- 



TABLE L Calculated lifetime broadening hr^^ in meV, of 
the n = 1 image state on Ag(lOO) with either ed{uj) — 1 
(or, equivalently, Zd — > —00) or £d{uj) from Ref. M] and 
Zd ~ — 1.5ao. Contributions to the linewidth from excita- 
tions where < uj < 3.5 eV are displayed in parentheses. The 
experimental linewidth is taken from time-resolved 2PPE ex- 
periments The lifetime in fs (Ifs — 10~^^ s) is obtained 
by noting that h — 658meVfs. 



Zd 


local part 


non-local part 


total 


experiment 




34(25) 


-16(-10) 


18(15) 




-1.5 


59(25) 


-47(-15) 


12(10) 


12 



citation channel which yields an overall enhancement of 
Im[-E]. 

Within a local picture of the self-energy and assuming, 
therefore, that Im[— E] near z w z' dominates the spec- 
trum, it follows from Eq. (|l|) that an overall enhancement 
of the self-energy yields an accordingly larger expectation 
value of Im[— E], i.e., an enhanced lifetime broadening. 
The key point of our work is, however, that this trend 
is reversed, due to the characteristic non-locality of the 
self-energy near the surface. As indicated in Table I, the 
calculated lifetime broadening is 12 meV, while without 
sd screening and without plasmon decay the broadening 
is much larger, namely, 18meV. 

To understand this puzzling result we note that ac- 
cording to Fig. 1 as long as z' < the self-energy exhibits 
local character. For z' > 0, however, non-local behavior 
dominates. To distinguish these two spatial regions, we 
refer to the broadening contribution arising from z, z' < 
as the 'local' part and the remainder {z > or z' > 0) 
as the 'non-local' part. These contributions to the life- 
time broadening are also given in Table I. We point out 
that non-local effects give rise to interference contribu- 
tions to the lifetime broadening, which are determined 
by the specific shape of the image-state wave function in 
the surface region where Im[— E] is maximum. 

This can be seen more clearly in Fig. ^ which shows the 
imaginary part of our calculated n = 1 image-state self- 
energy on Ag(lOO) (see also the bottom panel of Fig. |^), 
together with the n = 1 image-state wave function. Since 
this function has a node at z 1 oq the major inter- 
ference contributions to the lifetime broadening coming 
from z' > and z ^ z' must be negative. Further- 
more, as interference dominates the image-state relax- 
ation when z > or z' > 0, the lifetime broadening is 
actually smaller than one would have expected from the 
penetration of the image state. According to the results 
given in Table I, the role that d electrons play enhanc- 
ing the magnitude of the image-state self-energy is more 
pronounced on the vacuum side of the surface (where 
negative interference dominates) than in the bulk. The 
net effect of decay via surface plasmons and non-locality 
of the self-energy is therefore a considerably reduced life- 
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FIG. 2: Imaginary part of the n = 1 image-state self-energy 
Im[— Ek|| ,e{z, z')], versus z, in the vicinity of the (100) surface 
of Ag, for z' — 10 ao and Zd ~ — 1.5ao. The dotted line 
represents the n — 1 image-state wave function. As in Fig. |ll 
full circles represent the atomic positions and the value of z 
is indicated by an open circle, kii — and e = Ef + 3.9 eV. 



time broadening, in excellent agreement with the time- 
resolved 2PPE measurement reported in Ref . |Q . 

To investigate in more detail the impact of surface- 
plasmon excitation on the image-state lifetime broaden- 
ing, we also give in Table | the calculated linewidths ob- 
tained only from low-energy excitations {0 < u; < 3.5 eV) 
where the surface-plasmon relaxation channel is excluded 
(numbers in parentheses). In this case, the impact of 
d electrons on the lifetime broadening is much smaller. 
Thus, the opening of the surface-plasmon relaxation 
channel, which only occurs when both d electrons and all 
contributions with < w < 3.9 eV are included, yields 
a considerable enhancement of both the 'local' and the 
'non-local' part of the linewidth. As the large enhance- 
ment of the negative 'non-local' linewidth outweighs the 
increase of the 'local' contribution, the overall influence of 
the occupied d bands results in a reduced lifetime broad- 
ening. 

In conclusion, we have investigated the influence of 
single-particle and collective surface excitations on the 
decay of image-potential states on silver surfaces, by com- 
bining an accurate description of the dynamics of sp 
valence electrons with a physically motivated model in 
which the occupied d bands are accounted for by the 
presence of a polarizable medium. Our results demon- 
strate that despite the enhancement of the image-state 
self-energy due to decay via surface plasmons, the highly 
nonlocal character of this self-energy in the surface region 



ultimately leads to a surprisingly small lifetime broaden- 
ing, in agreement with the experimental data. 
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